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We present a parent Hamiltonian for weakly dimerized valence bond solid states for arbitrary half- 
integral S. While the model reduces for S = 1/2 to the Majumdar-Ghosh Hamiltonian we discuss 
this model and its properties for S = 3/2. Its degenerate ground state is the most popular toy 
model state for discussing dimerization in spin 3/2 chains. In particular, it describes the impurity 
induced dimer phase in Cr8Ni as proposed recently. We point out that the explicit construction of 
the Hamiltonian and its main features apply to arbitrary half-integral spin S. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Dg 



I. INTRODUCTION 



Valence bond solid (VBS) states played an impor- 
tant role for the understanding of quantum spin chains. 
Originally introduced as exact ground states of the inte- 
gral spin S Affleck-Kennedy-Lieb-Tasaki (AKLT) mod- 
els [HHI , VBS states meanwhile refer to all spin S states 
consisting of nearest-neighbor singlet bonds. Amongst 
them, the spin 1/2 Majumdar-Ghosh (MG) states 0,11] 
are the most important examples as their parent Hamil- 
tonian corresponds to a certain point in the parameter 
space of the spin 1/2 J\-Ji model with antiferromagnetic 
Jj. As such, the MG model links dimerized VBS states 
and dimerized spin chains where the frustration is forced 
by next-nearest neighbor interactions. For S > 1/2, 
however, this exact link is not given anymore. The VBS 
states seem not to belong to the phase diagrams of the 
frustrated Heisenberg chains. Nonetheless, VBS states 
are often used as toy model states in order to discuss 
frustration in quantum spin chains, see for the spin 3/2 
case for example Refs. @, H 0, 1, M, EH, GJ, [ll and 
references therein. Parent Hamiltonians, however, are 
only known for the Majumdar-Ghosh states for the 
translational invariant AKLT states [l[ , and for the spin 
1 dimer states p|. For spin 3/2 and larger S, no exact 
and translational invariant Hamiltonians are known hav- 
ing dimerized VBS states as their unique ground states. 
Most recently, VBS states have attracted renewed in- 
terest. They have been generalized to special unitary 
SU(n) symmetry 



symmetr 
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161 . 1171 ] , to special orthogonal 
19l | . and to symplectic symmetry 
A supersymmetric extension of VBS states 
has been proposed recently pH ]. 

In this Letter, we present a parent Hamiltonian for the 
partially or weakly dimerized VBS states (see Fig.Q] for 
the illustration of the spin 3/2 case). For spin 1/2, this 
Hamiltonian reduces to the MG Hamiltonian. The par- 
ticular importance of the spin 3/2 Hamiltonian is given 
by its possible realization in the impurity induced dimer 
phase in Cr8Ni as proposed recently [12| ■ This molecu- 
lar nanomagnet comprises eight chromium(III) ions with 
spin 3/2 and one nickel(II) ion with spin 1. The ground 
state of this Cr8Ni ring was shown to be well understood 
in the VBS picture [12| by means of numerical Simula- 



FIG. 1: Illustration of the spin 3/2 dimer states, (a) and (b) 
are the two-fold degenerate partially dimerized VBS states. 
Each state can be seen as the covering of the AKLT state with 
one of the Majumdar-Ghosh states. The large grey-shaded 
circles are the lattice sites, while the small black dots denote 
a spin 1/2. The spins 1/2 on a lattice site will be projected 
onto the totally symmetric subspace, i.e., S = 3/2. The lines 
between two spins 1/2 denote a singlet bond. 



tions. Likewise, the Hamiltonian can be implemented 
within DMRG, as only standard spin interaction terms 
are present. This will allow for extensive numerical stud- 
ies complementary to known results for the spin 3/2 J\- 
Ji model with antiferromagnetic Ji @. 

The Letter is organized as follows: we start with a brief 
discussion of general properties of VBS states. Then we 
introduce the S — 3/2 dimer Hamiltonian which has the 
weakly dimerized VBS states shown in Fig.Q] as ground 
states and show that these are the unique ground states. 
We show that the lowest lying excitations are fraction- 
ally quantized spin 1/2 spinons. Then we discuss the 
connection to the dimerized phase in the frustrated spin 
3/2 Heisenberg model and static spin-spin correlations. 
Finally, we show that the explicit construction of the 
Hamiltonian as well as all qualitative results generalize 
to the case of arbitrary half-integral spin S. 



II. VBS STATES 

In 1987, Affleck, Kennedy, Lieb and Tasaki (AKLT) 
introduced a spin 1 Hamiltonian which was constructed 
as a sum over a local projection operator. Its unique 
ground state can be understood as follows: on each lat- 
tice site are two spins 1/2 symmetrically coupled into an 
effective spin 1, while each of the spins 1/2 is antisym- 
metrically coupled to a spin 1/2 on its neighboring lattice 
sites, either to the right or the left. For an illustration see 
Fig.[U(a). This state was called a valence bond solid. We 
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FIG. 2: The 5 = 1 VBS states: (a) is the translational in- 
variant AKLT state and (b) is one the of two-fold degenerate 
dimer states. The second dimer state can be obtained by 
shifting the state (b) by one lattice spacing to the right. 



wish to understand as a VBS all states consisting of local 
singlet bonds. In particular, we call the MG states and 
the spin 1 dimer states (see Fig. [5] (b)) also VBS states 
even though they break translational invariance sponta- 
neously. 

One of the main features of VBS states is that we are 
able to write down their wave functions explicitly as eas- 
ily as we illustrate them. Schwinger bosons [13, [H, [24| 
constitute a convenient way to formulate spin 5 repre- 
sentations. Arovas et al. [25[ applied for the first time 
Schwinger bosons to the AKLT model and to VBS states 
in general. The spin 5 operators 



5+ = a^b, S~ = tfa, 



a) a 



btft) (1) 



are expressed by bosonic creation and annihilation op- 
erators which obey the standard commutation relations 
[b, f>T] = l while all other commutators vanish. 



[a, a 



tl - 



The spin quantum number S is given by 25 = a^a + b^b 
and the general spin 5 state is defined as 



FIG. 3: Illustration of one of the two-fold degenerate com- 
pletely dimerized S = 3/2 VBS states, say \ip(o,3))- By shift- 
ing of one lattice spacing to the right we obtain the second 
dimer state |t/>(3,o))- 



(V'fSjS))' Note that there is always a pair of completely 
dimerized VBS states regardless 5. We stress again that 
the Schwinger bosons ensures the symmetric projection 
on every lattice site. 

In this Letter, we are mainly interested in 5 = 3/2 
VBS states. There arc four of them: completely dimer- 
ized states | "0(3,0) ) an d | "0(0,3)) (the latter is shown in 
Fig. [3]) and weakly dimerized states 10(2,1)) and |V>(i,2)) 
(both are shown in Fig.[T]). For the weakly dimerized 
VBS states we present in the next paragraph a parent 
Hamiltonian. We wish to mention that one can easily 
define a spin-5 Hamiltonian in terms of local projection 
operators annihilating all 25 + 1 VBS states. It involves 
three-site interactions and reduces for S = 1/2 to the 
well-known MG Hamiltonian. For S > 1/2, however, 
this is not a parent Hamiltonian anymore as it has all 
spin 5 VBS states as groundstates. The Hamiltonian is 
given by 



(4) 
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The spin 1 /2 states are thus given by cU | ) = a\ | ) and 

cjjj ) = b\ | ) , respectively. The difference between the 
fermionic operators and Schwinger bosons shows up only 
when two or more creation operators act on the same site. 
While fermion operators create a singlet configuration 
c 1t c Ii Schwinger bosons create a totally symmetric 
representation. 

The elegance of Schwinger bosons come to light when 
defining the general spin S VBS state: 



(3) 



where Bi = a\b\ +1 — b\a\ +1 creates on sites i and i + 1 a 
singlet bond. We recover for (m = l,n = 0) and (0,1), 
respectively, the MG states. The 5=1 AKLT state can 
thus be written as |V>(i,i)) and the spin 1 dimer state 
(Fig-Hi (b)) as |0(2.o))- Obviously, m + n = 2S and there 
are always 25+1 VBS states. For half-integer 5, there 
are always 5 + \ pairs of dimerized states with different 
strengths of dimcrization. For integer 5, there are 5 
pairs of dimerized states and one additional AKLT state 



III. SPIN 3/2 DIMER HAMILTONIAN 

Now we present the Hamiltonian which uniquely sin- 
gles out the weakly dimerized VBS states |'0(2.i)) and 
1 0(1,2)) as shown in Fig.[TJ This Hamiltonian is a sum of 

(S) 

projection operators i^ t+i i+ti wn ich project onto the 
spin 5 subspace on /i + 1 neighboring sites. 



n - V f M+M+2 + A i r i,i+1 



i+2 



Remarkably, the last projection operator acts on two 
neighboring sites while the other operators act on three 
neighboring sites. We decompose, hence, the Hamilto- 
nian into two parts, Hi acting on three neighboring sites 
and H.2 acting on neighboring sites. We are able to ex- 
press Tii = Wi(|) by means of spin operators as follows, 



n 1 = J2[(s l + s i+1 + s. l+2 ) 2 -- 



(5, + s l+1 + s. l+2 f - H 



(5) 



2 



The second part, H2, corresponds to the nearest neigh- 
bor bilinear-biquadratic-bicubic Hciscnbcrg model with 
certain parameters, 



Hi 



E 



Vp (3) 



(6) 



Finally, the parent Hamiltonian can be written for /3 > 
as 



(7) 



Note that this Hamiltonian can be implemented both 
within exact diagonalization (ED) and DMRG. Terms 
like (SiSj) 2 and (SiSj) 3 are standard terms. More com- 
plicated terms like (SiSj)(SiSk) or (SiSj)(SjSk) are 
also straight forward to implement. The reader may no- 
tice that DMRG has been successfully applied by several 
authors [H, [2?], [H, H^, [3(| to investigate antiferromag- 
netic spin 3/2 chains. 

In the following we will prove that ([7]) has a two-fold 
degenerate zero-energy ground state, 



H 



Ka.1)) =H dimoi > (1 , 2) ) = 0. 



(8) 



The proof is exceedingly simple: first we will verify that 
Hi has the four 5 = 3/2 VBS states as unique zero- 
energy groundstates, second we will show that H2 ex- 
hibits a huge (actually an infinite) ground state degen- 
eracy (also with zero energy) which contains the weakly 
dimerized VBS states but not the completely dimerized 
VBS states. Consequently, the additive term /3H.2 in ([7]) 
does not affect the weakly dimerized VBS state while the 
fully dimerized state will be lifted to higher energy. This 
lifting is controlled by the parameter (3 > and can be 
tuned continuously. 

The auxiliary Hamiltonian Tti has the four S = 
3/2 VBS states as exact zero-energy ground states, 

Hi 1^(3.0)) = Hi 1^(2,1)) = Hi I "0(1,2) ) = Hi |^(0,3)) = 

0. 

This can be seen as follows: Consider three neighbor- 
ing sites of one of the four VBS states (no matter which of 
the VBS states, completely or partially dimerized ones): 
there are always three singlet bonds and three individual 
spin 1/2 representations, 0(g)0(g>0(gi|® \®\ = 2-|ffi§. 
In general, the tensor decomposition of three spin 3/2 is 
givenby§<g>§g>§ = 2-|©4-§ffi3-§ffi2-§ffi§. The pro- 
jection onto the subspace with representation S = 5/2, 
S = 7/2, and S = 9/2 singles out the four VBS states 
as ground states. We can express the local projection 
operator in terms of spin operators as done in equation 
(|5"|). Now we consider the nearest neighbor projection 



(3) 

operator Pl i+1 - The Hamiltonian H2 = Y^i + a+i 
an infinite ground state degeneracy. This can be seen as 
every state is a ground state where two spins 1/2 per site 
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FIG. 4: Lowest lying energies of the Hamiltonian (ffj) as a 
function of /3 at the T-point (momentum k = 0) in the Bril- 
louin zone. The data was obtained on a chain with N = 12 
lattice sites and periodic boundary conditions. As the T-point 
is considered, there are only two of the four groundstates for 
/3 = 0. The other groundstates ( |V>(i,2)) and |V'(o.3)) ) are lo- 
cated at k = it /a in the Brillouin zone. For increasing f3 > 
the fully dimerized VBS state j "0(3,0)) is increased in energy. 
Around f3 ~ 1.1 it crosses with the first excited state of Tti 
which is in some sense also the first excited state of 7-^ dlmer . 
This crossing is expected, but the value f3 ~ 1.1 depends on 
the system size and will decrease for larger N. 



are in the "AKLT configuration" and the third spin 1/2 
is "free" . The partially dimerized states obviously belong 
to this set of ground states while the completely dimer- 
ized states are not ground states of H2- Consequently, 
the term f3H2 in (J7]) lifts the completely dimerized states 
continuously to higher energy as /3 is increased. The par- 
tially dimerized states remain as unique ground states. 
This scenario is shown in Fig.[J]for a chain with N = 12 
sites and periodic boundary conditions where the data 
was obtained by means of exact diagonalization. 

For open boundary conditions the situation changes: 
the dimer states are no global singlets anymore, edge 
states emerge, and the dimer states arc fourfold (sin- 
glet and triplet) and ninefold (singlet, triplet, and quin- 
tuplet) degenerate. When cutting the dimer states, ei- 
ther one individual spin 1 /2 or two symmetrically cou- 
pled spins 1/2 remain which are not coupled in local sin- 
glet bonds (the two different cuts correspond to the two 
ground states in case of periodic boundary conditions). 
The dangling spins at the edges cause the degeneracy of 
the dimer ground states. For the first situation (A^ even) 
the two edge spins couple into a singlet and a triplet, 
\ <8> \ = © 1, yielding a four-fold degenerate dimer 
state. For the second situation (N even), the two spins 
1/2 at each edge are symmetrically coupled into a triplet, 
5{g ® |} = 1- Hence, the two edge spins carry spin 1 
and will be coupled into singlet, triplet, and quintuplet, 
101 = 0©1©2, yielding a ninefold degenerate ground 
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FIG. 5: (a) The only possibility to excite the system being in 
one of the VBS ground states is to break a singlet bond (indi- 
cated by the dashed line), (b) The resulting spinons (the small 
dark-grey dots) are deconfined and can freely move through 
the chain as they interpolate between different ground states. 
In fact, the left spinon in the cartoon constitutes a domain 
wall between the ground states V(2,i)) ant ^ | "0(3,0) ) while the 
right spinon constitutes a domain wall between |V>(3,o)) an d 
(2,i))- ( c ) For clarity we have fixed the left spinon. From 
the sequence (a), (b), and (c) we note that the spinons can 
occupy only half the number of sites (see the explanation in 
the text). The reader should notice that the "states" shown 
in (a), (b), and (c) are not exact eigenstates of Tii but the 
cartoons provide a good intuitive picture of the low-lying ex- 
citations. 

state. This manifold degeneracy of the ground state in 
case of open boundary conditions can be nicely observed 
numerically, when looking at the different S^-subspaces. 



FIG. 6: The spinon excitation (the small dark-grey dot) in- 
terpolates between the states |V'(2,i)) an d 1^(1,2) )• The state 
on the right hand side of the spinon excitations is identical to 
the state in Fig.[T]b as explained in the text. 

is irrelevant. Hence, a VBS state is identified uniquely by 
to and n. The cartoon in Fig. [5] suggests that elementary 
excitations of {7J) are also spin 1/2 spinons. This is not 
unexpected as the Hamiltonian ([JJ is the Hamiltonian Hi 
plus an additional term 7i 2 which lifts up the fully dimcr- 
ized states out of the ground state manifold. Altogether, 
the whole consideration of deconfined spinon excitations 
in the spin 3/2 dimer model is consistent with the picture 
that the origin of the Haldane gap [HJ is a confinement 
force between spinons [T^, [3^] . The generic S = 3/2 spin 
chain, the nearest neighbor Heiscnbcrg model, is known 
not to exhibit a Haldane gap and, hence, the spinons 
cannot be confined. While the generic spin 3/2 chain has 
a gapless excitation spectrum, the Hamiltonians 7-^ dlmcr 
and Hi as well as all spin S Hamiltonians mentioned at 
the beginning of this section exhibit a gapped excitation 
spectrum as it costs a finite energy to break a singlet 
bond. This energy cost survives in the thermodynamic 
limit but it is not Haldane's gap. 



IV. FRACTIONALLY QUANTIZED 
EXCITATIONS 



As already discussed, we are able to define the model 
Hamiltonian Tti(S) (see Eq. Qj) acting on three neigh- 
boring sites for arbitrary spin S exhibiting a (2 5+1) -fold 
degenerate ground state, i.e., the 25 + 1 possible spin 5 
VBS states |hJ S3 • For 5=1/2, this model corresponds 
to the MG Hamiltonian, for S = 1 to the Hamiltonian 
proposed by Sen and Surendran 31] , and for 5 — 3/2 to 



Hi of equation In all of these Hamiltonians, the low- 
est lying excitation will be a pair of deconfined 5=1/2 
spinons. This is easy to understand as the spinons emerge 
in this models as domain walls interpolating between dif- 
ferent ground states. As the Hilbert space of spin 5 an- 
tifcrromagnets is spanned by spin flips which carry spin 
1, the spin 1/2 excitation has to be seen as a fractionally 
quantized excitation. Wc illustrate this situation for the 
Hamiltonian Hi in Fig. [5] In the cartoon Fig.[5l we no- 
tice another interesting point: even though the spinons 
are free they can only occupy every second lattice site. 
If there are only a few spinons in a chain the number of 
orbitals available to them is roughly half the number of 
sites reflecting their fractionally quantized character. 

Now the question remains whether this picture of free 
spinon excitations holds for the Hamiltonian 7^ dlmcr con- 
sidered in this paper. As we mentioned above, relevant 
in our cartoons are only the number of singlet bonds be- 
tween neighboring sites while the order of singlet bonds 



V. DISCUSSION 

The Hamiltonian ([7]) has the weakly dimerized states 
1^(2,1)) and 2))i respectively, as exact ground states. 
Dimerization is usually caused by frustration. In one spa- 
tial dimension this is realized in a next-nearest neighbor 
Heisenberg model, H = J2i + a &iSi+2, when the 

next nearest neighbor coupling a exceeds a critical value 
etc- 

The spin 1/2 Heisenberg model is characterized by a 
critical phase for a < a c = 0.2411 [H [H, [H. The 
system undergoes a second-order phase transition to a 
dimerized phase for a > a c . For cvmg, the Hamiltonian 
corresponds to the Majumdar-Ghosh model exhibiting 
the exact dimer states k/vi.o)) and jV^o,!)) as ground 
states. Notice, that the dimerization becomes exact in 
the thermodynamic limit only. That is, the ground state 
becomes perfectly degenerate due to the broken transla- 
tion symmetry. One ground state has total momentum 
k = while the other has k = ir. Except for cvmg ■ the 
ground state degeneracy is not intact in the finite sys- 
tem. Nevertheless, ground states in the dimer phase are 
known to exhibit a good overlap with the MG ground 
states. 

For the spin 3/2 Heisenberg model with next near- 
est neighbor interaction a, Roth and Schollwock found 
the phase transition to the dimerized phase at a c ^/ 2 ~ 
0.29 0]. They further computed the maximum dimeriza- 
tion for cv max = 0.415. We expect that the ground states 



4 



in the spin 3/2 dimcr phase correspond to the weakly 
dimerized states rather than the completely dimcrized 
states, i.e., our Hamiltonian ([7]) can be seen as the 
S = 3/2-analogue of the MG model. We can fur- 
ther conjecture that the states V(2,i)) an d ^(l^)) wm 
have an excellent overlap with the ground states in the 
dimer phase investigated by Roth and Schollwock for 
sufficiently large system lengths. Such system lengths 
are not computable within ED. Nevertheless, in princi- 
ple the Hamiltonian ([7| is implcmentable within DMRG 
as pointed out above and the calculation of overlaps 
might be feasible. Furthermore, it might be interest- 
ing to compute the static spin-spin correlations within 
DMRG. Correlations are known to decay quite fast in 
dimer phases, for the MG states, however, the correla- 
tions decay abruptly, i.e., {SiSi+2} = 0. Hence, the 
question remains for the spin 3/2 dimer model whether 
the spin-spin correlations decay fast, say exponentially, 
or abruptly. 



dimerized" states coincide with the completely dimcr- 
ized states). Since the second sum in ^ does not con- 
tribute for S = 1/2 the Hamiltonian (J9j> reduces to the 
sum over the projection operator onto the S = 3/2 sub- 
space on three neighboring sites which is known to be 
the MG Hamiltonian. For S = 3/2, the Hamiltonian 
TL S is equivalent to ([7|). For larger S, e.g. S = 5/2, 
we find easily the Hamiltonian Tii = which has 

a 2S + 1-fold degenerate ground state manifold consist- 
ing of all VBS states. Again the term T-i 2 lifts up all 
VBS to higher energy except the weakly dimcrized states. 
For S = 5/2, the Hamiltonian H2 is given by H2 = 
£J 13397(5^+1) + 3582(S 4 S m ) 2 + 400(S l S 4+1 ) 3 + 
16(SiSi +1 ) 4 + 274505/16]. Although one is restricted 
to small chain lengths when S becomes 5/2 or larger, 
all these terms are implementable within ED and the 
5 = 5/2 spin-dimer model Hi + PH2 might be verified 
numerically. 

VII. CONCLUSION 



VI. GENERALIZATION TO ARBITRARY 
HALF-INTEGER S 

We are able to formulate the generalized version of the 
Hamiltonian j7| to arbitrary half integer S. It involves 
interactions between three neighboring sites and is for- 
mulated in terms of projection operators as follows: 



H 




A P (m) 

-^Mi r i,i+l,i+2 
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2=S+ 



A P (/X2) 



(9) 

We have restricted S to be half-integer and (3 > 0. 
The two-fold degenerate ground state of H s is explicitly 

given by "0(s±±,s=f±)^ by means of ([3]). For S = 1/2 
the weakly dimerized states reduce to the well-known 
Majumdar-Ghosh (MG) states (in this case, the "weakly 



In conclusion, we have presented a parent Hamilto- 
nian for the weakly dimerized VBS states with arbitrary 
half-integral 5. For spin 3/2, we have explained its con- 
struction explicitly and have discussed its main proper- 
ties. These weakly dimerized spin 3/2 VBS states might 
be realized in the impurity induced dimer phase of the 
ground state of the molecular nanomagnct Cr8Ni. The 
path for further numerical consideration has been pointed 
out. Finally, the generalization to spin S has been ex- 
plained. 



VIII. ACKNOWLEDGEMENT 

I am grateful to Martin Greiter for numerous discus- 
sions and his encouragement for this paper. I am also 
indebted to Peter Wolfle and Hong-Hao Tu for several 
helpful suggestions. 



[1] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. 

Rev. Lett. 59, 799 (1987). 
[2] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Com- 

mun. Math. Phys. 115, 477 (1988). 
[3] C. K. Majumdar and D. K. Ghosh, J. Math. Phys. 10, 

1388 (1969). 

[4] C. K. Majumdar and D. K. Ghosh, J. Math. Phys. 10, 
1399 (1969). 

[5] I. Affleck, J. Phys.: Condens. Matter 1, 3047 (1989). 
[6] D. Guo, T. Kennedy, and S. Mazumdar, Phys. Rev. B 

41, R9592 (1990). 
[7] H. Niggemann and J. Zittartz, Z. Phys. B 101, 289 

(1996). 

[8] M. Oshikawa, M. Yamanaka, and I. Affleck, Phys. Rev. 

Lett. 78, 1984 (1997). 
[9] R. Roth and U. Schollwock, Phys. Rev. B 58, 9264 

(1998). 



[10] M. Nakamura and S. Todo, Phys. Rev. Lett. 89, 077204 
(2002). 

[11] J. Almeida, M. A. Martin-Delgado, and G. Sierra, Phys. 

Rev. B 77, 094415 (2008). 
[12] J. Almeida, M. A. Martin-Delgado, and G. Sierra, Phys. 

Rev. B 79, 115141 (2009). 
[13] G. Refael, S. Kehrein, and D. S. Fisher, Phys. Rev. B 66, 

060402 (R) (2002). 
[14] S. Chen, C. Wu, S.-C. Zhang, and Y. Wang, Phys. Rev. 

B 72, 214428 (2005). 
[15] M. Greiter, S. Rachel, and D. Schuricht, Phys. Rev. B 

75, 060401 (R) (2007). 
[16] M. Greiter and S. Rachel, Phys. Rev. B 75, 188841 

(2007). 

[17] D. P. Arovas, Phys. Rev. B 77, 104404 (2008). 
[18] H.-H. Tu and G.-M. Zhang and T. Xiang, J. Phys. A: 
Math. Theor. 41 (2008) 415201 



5 



[19] H.-H. Tu and G.-M. Zhang and T. Xiang, Phys. Rev. B 

78, 094404 (2008). 
[20] D. Schuricht and S. Rachel, Phys. Rev. B 78, 014430 

(2008). 

[21] D. P. A rovas, K. Hasebe , X.-L. Qi, and S.-C. Zhang, 

preprint larXiv:0901.1498l (2009) . 
[22] J. Schwinger in L. Biedenharn and H. van Dam (eds), 

quantum Theory of Angular Momentum (Academic 

Press, 1965). 

[23] A. Auerbach, Interacting electrons and quantum mag- 
netism (Springer, New York, 1994). 

[24] D. P. Arovas and A. Auerbach, preprint larXiv:0809.4836l 
(2008). 

[25] D. P. Arovas, A. Auerbach, and F. D. M. Haldane, Phys. 

Rev. Lett. 60, 531 (1988). 
[26] S. Qin, T.-K. Ng, Z.-B. Su, Phys. Rev. B 52, 12844 

(1995). 

[27] K. Hallberg, X. Q. Wang, P. Horsch, and A. Moreo, Phys. 



Rev. Lett. 76, 4955 (1996). 
[28] E. Carlon, P. Lajko, H. Rieger, F. Igloi, Phys. Rev. B 69, 
144416 (2004). 

[29] G. Fath, 6. Legeza, P. Lajko, and F. Igloi, Phys. Rev. B 

73, 214447 (2006). 
[30] F. Heidrich-Meisner, I. A. Sergienko, A. E. Feiguin, and 

E. R. Dagotto, Phys. Rev. B 75, 064413 (2007). 
[31] D. Sen and N. Surendran, Phys. Rev. B 75, 104411 

(2007). 

[32] F. D. M. Haldane, Phys. Lett. A 93, 464 (1983); Phys. 

Rev. Lett. 50, 1153 (1983). 
[33] M. Greiter, J. Low Temp. Phys. 126, 1029 (2002). 
[34] K. Okamoto and K. Nomura, Phys. Lett. A 169, 433 

(1992). 

[35] S. Eggert, Phys. Rev. B 54, R9612 (1996) 
[36] R. Jullien and F. D. M. Haldane, Bull. Am. Phys. Soc. 
28, 344 (1983). 



G 



